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We have found a novel feature of the bistable large-spin model described by the Hamiltonian Ti = 
—DS^ — HxSx- The crossover from thermal to quantum regime for the escape rate can be either first 
{Hx < SD/2) or second (SD/2 < Hx < 2SD) order, that is, sharp or smooth, depending on the 
strength of the transverse field. This prediction can be tested experimentally in molecular magnets 
like Mni2Ac. 

PACS numbers: 75.45. -fj, 75.50.Tt 



Transitions between two states in a bistable system 
can occur either due to the classical thermal activation 
or via quantum tunneling. A rigorous study of that prob- 
lem was begun by Kramers and WKB and a 
review of the progress that followed can be found in Ref. 
1^. At high temperature the transition rate follows the 
Arrhenius law, T ~ exp(— A[//T), with AU being the 
hight of the energy barrier between the two states. In 
the limit of T ^ 0, the transitions are purely quantum, 
r ^ exp(— _B), with B independent on temperature. Due 
to the exponential dependence of the thermal rate on T, 
the temperature Tq of the crossover from quantum to 
thermal regime can be estimated as T^°^ = AU/B. For 
a quasiclassical particle in a potential U{x), Goldanskii 
1^ noticed the possibility of a more accurate definition, 
Tq — h/TQ, where tq is the period of small oscillations 
near the bottom of the inverted potential, ~U{x). Be- 

low Tq , thermally assisted tunneling occurs from the 
excited levels, that reduces to the tunneling from the 
ground-state level at T = 0. Above Tg ' quantum effects 
are small and the transitions occur due to the thermal 
activation to the top of the barrier. Affleck 0| demon- 
strated that the two regimes regimes smoothly join at 
T = Tq^-*. Larkin and Ovchinnikov called this sit- 
uation the second-order phase transition from classical 
to quantum behavior. This means that for F written as 
F ~ exp(— AC//Toff ), the dependence of both Toff and its 

(2) 

first derivative on T are continuous at T = Tg . This 
situation is not generic, however. The transition between 
the two regimes can also be of the first order i.e., 
more abrupt, with dTes/dT discontinuous at a certain 
temperature Tq^-* > T^p . Chudnovsky derived the cri- 
terium allowing one to establish whether first- or second- 
order transition takes place, based on the shape of the po- 
tential U{x). Commonly studied potentials U = —x'^+x* 



and U 



x^ yield the second-order transition. 



Physically relevant potentials which would exhibit the 
first-order transitions were not known. In this work we 
show that spin systems readily accessible in the experi- 
ment possess both first- and second-order transitions be- 
tween the classical and quantum behavior of the escape 
rate. The order of the transition in these systems can be 
controlled by external magnetic field. 

Consider a spin system described by the Hamiltonian 



(1) 



where S ^ 1. This model is generic for prob- 
lems of spin tunneling studied by different methods 
0,0jl2[|l|,|l|,|l|l. It is believed to be a good approxi- 
mation for the molecular magnet Mni2 Ac of spin S" = 10, 
intensively studied in last years (see, e.g., Ref. |jl^). In 
the quasiclassical approximation the transition rate is 
given by 



dE I^(i?)e-(^-^-")/^, (2) 

where W{E) is the probability of tunneling at an energy 
E and i?min corresponds to the bottom of the potential. 
This probability is defined via the imaginary-time action 

W[E) - e-^(^). (3) 
With the accuracy to the exponent, 

F^e-WT^ (4) 

where i^min is the minimum of the effective "free energy" 

F=E + TS{E) - Sniin (5) 

with respect to E. 

In order to obtain S{E) for the Hamiltonian (||) we 
will use the method of mapping the spin problem onto 



1 



a particle problem ||lj,|lj,y/7|,^ . The equivalent particle 
Hamiltonian is 



V2 

n^~— + uix), 

2m 



(6) 



where 



and 



- 1 _ 

m = , Air = ^ ^ ■ 

2D {2S+1)D 



(8) 



In the future we shall neglect 1/2 in comparison to 5 ^ 1. 

The imaginary-time action is then given by the WKB 
expression 



x{E) 



S{E) = 2(2m)i/2 J dx ^U{x) - E, 



(9) 



-x(E) 



where ±x{E) are the turning points for the particle oscil- 
lating inside the inverted potential —U{x). The period 
of these oscilations, Tp{E) = —dS{E)/dE, depends on 
energy. Minimization of (o) gives 



(10) 



the condition familiar from the quantum statistics 
p9|,[7|P] . It determines the instanton trajectory that dom- 
inates the transition rate at a temperature T. 

The dependence of Tp{E) on E determines the kind of 
the crossover from quantum tunneling to thermal activa- 
tion . If Tp monotonically increases with the amplitude 
of oscillations, i.e., with decreasing energy E, the tran- 
sition is of the second order. This kind of the crossover 
has been intensively studied, including the case of tun- 
neling with dissipation |2^,|2^,^,^,|^ . If, however, the 
dependence of Tp{E) is non-monotonic, the first-order 
crossover takes place. Let us demonstrate that both 
kinds of the crossover exist for our spin model, depend- 
ing on the strength of the transverse field. Expanding 
near x = 0, one obtains 



U{x) ^ U{0) + S'^D 

+ — 
3 



hx - ^ \ X 



hx{l - K)x^ 
' + 0{x^) 



(11) 



where the sixth-order term is positive. The second-order 
term in ( p]) is negative for hx < 1, which corresponds to 
the existence of the energy barrier 



Umax t^min — S^D{1 hx)"^ 



(12) 



For hx > 1/4 the fourth-order term in ( [ll| ) is positive. 



i.e., U{x) is of the form —x 



The inverted potential 



— U{x) is hence of the type x'^ — x*, which results in the 
increase of Tp with the oscillation amplitude (i.e., with 
lowering the energy E) and to the second-order transi- 
tion Q. At /la; < 1/4 the anharmonicity of —U{x) has 
the opposite sign, —U{x)'-^x^ + x'^, which leads to the 
decrease of Tp when lowering E for energies below the 
top of the barrier. However, with further lowering of E 
the period Tp begins to increase and diverges logarith- 
mically for E approaching the bottom of the potential. 
This non-monotonic behavior of Tp [E) leads to the first- 
order transition from the thermally activated escape to 
the quantum escape |^ . 

In the case of the second-order transition the crossover 
occurs at temperature 



(2) _ t^O 



SD 



TT- = \/hx{l - hx 

2n TT 



(13) 



where Cjq ^ ^/\U"{0)\/ m is the instanton frequency 
It is convenient to introduce dimensionless temperature 
and energy variables 



(2)^ 



P = 



(14) 



The effective free energy ^ near the top of the barrier 
(P << 1) can be calculated with the use of Eqs. (^) and 
( prj ) and reads 



FiP) 



^ 1 



1 



- 1)P + - ( 1 —] 

' 8 V 4/i, ' 



1 



(15) 



The analogy with the Landau theory of phase transitions, 
described hy F = atp'^ + hij]'^ -t- , now becomes appar- 
ent. The factor in front of P (the Landau coefficient 
a) changes the sign at the phase transition temperature 
T = T^^\ The factor in front of (the Landau coef- 
ficient b) changes the sign at the field value hx = 1/4 
determining the phase boundary between the first- and 
the second-order transitions, as has been already noticed 
from Eq. (pi]). The factor in front of P^ (the Landau 
coefficient c) remains always positive. The numerically 
computed dependence of P on P for the entire range of 
energy is plotted in Fig. pi. At hx = 0.3 (Fig. la) the 

• (2) 

minimum of F remains J/max ^ C^min for all T > Tq . Be- 

(2) 

low Tg ' it continuously shifts from the top to the bottom 
of the potential as temperature is lowered. This corre- 
sponds to the second-order transition from thermal acti- 
vation to thermally assisted tunneling. At hx — 0.1 (Fig. 
2b), however, there can be one or two minima of P, de- 
pending on temperature. The crossover between classical 
and quantum regimes occures when the two minima have 
the same free energy, which for hx = 0-1 takes place at 

The crossover temperature for the escape rate is fre- 
quently estimated by equating the ground-state tunnel- 
ing exponent to that of thermal activation: 
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FIG. 1. EfTective "free energy" for the escape rate: (a) 
- hx = Hx / {2SD) = 0.3, second-order transition; (b) - 
hx = 0.1, first-order transition. 
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The ground-state tunneling exponent B given by Eq. (|^) 
can be analytically calculated jl^, which together with 
Eq. (|l|) yields 



(0) 



SD 



(1 



4 Jii^JIEK 



SI^l ln[2/(e;i,)]' 



K < 1 



(17) 



One can see from Fig. lb that T^^ somewhat under- 
estimates the crossover temperature. For hx = 0.1 one 
has T^°^ = l.Oeird^^ < T^a^\ The estimation T^^^ be- 
comes, however, accurate in the limit of small hx- The 
dependence of the crossover temperature To on the trans- 
verse field in the whole range, < /ix < 1, is pre- 
sented in Fig. 2. The temperature dependence of the 
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FIG. 2. Dependence of the crossover temperature To on the 
transverse field. 




FIG. 3. Dependences of the effective temperature TcS on 
T for the different values of the transverse field. 



escape rate can be conveniently written in the form 
r ~ exp[— (f/inax — f^min)/7cff (T)], whcrc the dependence 
of the effective temperature on T is presented in Fig. 3 
for different hx ■ It can be seen from Fig. 3 that the most 
significant difference between the crossover temperature 
Tq"' of Eq. (^6|) and the actual crossover temperature Tq 
arises in the limit of small barrier, that is, at hx — *■ 1. 
The former is described by the intersection of the dotted 
Arrhenius line with the horizontal line corresponding to 
Toff (T)/To at zero temperature. From Eqs. (|l|) and (|l|) 
for hx 1 one obtains T^°^ /T^j'^^ = 3n/{8V2) « 0.833. 

As follows from Fig. 3, the difference between the 
curves Teff(T) describing the first- and second-order 
crossover is quite dramatic. It must be easily observed 
in experiment if the appropriate system is found. Very 
recently experiments on individual small magnetic par- 
ticles with S ^ 10^ — 10^ have become possible pS]. In 
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these experiments the barrier is lowered by tuning the 
magnetic field to the critical value. In our model this is 
the case of the second-order transition. In order to get 
the first-order transition, must be lower than Ha/ 4:, 
where Ha = 2SD/{g^B) is the anisotropy field. This 
case requires a moderate spin S in order to provide a 
significant escape rate. The Hamiltonian has been 
found to be a good model for Mni2Ac S — 10. In 
this case the quantization of spin levels becomes impor- 
tant. However, our statement regarding the possibility 
of first- and second-order transitions remains valid . 

The analogy with phase transitions in the temperature 
dependence of the escape rate formally exists only in the 
limit of S* ^ oo. For a finite 5", the transition from (||) to 
has the accuracy of 1/ S. Quantum corrections to the 
escape rate above Tq are of the same order. Thermal 
and quantum corrections will smoothen the first-order 
transition in the narrow remperature region close to Tq. 
Nevertheless, even for S = 10, the difference between the 
crossover at small and large H^ must be easily observ- 
able. The sharpness of the crossover between thermal 
and quantum regimes also depends on the strength of 
the dissipation. In the case of the low dissipation which 
is common for the magnetic systems, its effect on the 
crossover is smah §. 

For Mni2Ac the anisotropy field is about 10 T. The 
crossover from thermal to quantum regime should, there- 
fore, switch from first to second order at H^ ~ 2.5 T. The 
crossover temperature is about 1 K |^^. These ranges 
of field and temperature are easily accessible in exper- 
iment. Note that similar effects may exist in the Feg 
molecular magnet where the crossover from thermal to 
quantum regime has been already observed This 
system, however, is described by the spin Hamiltonian 
with the transverse anisotropy which requires separate 
theoretical investigation. We believe that the statement 
made in this paper, regarding the possibility of first- and 
second-order crossover from thermal to quantum regime, 
must be very general for spin systems. 
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